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ABSTRACT 

We study the orbits, tidal heating and mass loss from satellites around close-in gas giant exoplanets. 
The focus is on large satellites which are potentially observable by their transit signature. We argue 
that even Earth-size satellites around hot Jupiters may be immune to destruction by orbital decay; 
detection of such a massive satellite would strongly constrain theories of tidal dissipation in gas giants, 
in a manner complementary to orbital circularization. The star's gravity induces significant periodic 
eccentricity in the satellite's orbit. The resulting tidal heating rates, per unit mass, are far in excess of 
lo's and dominate radioactive heating out to planet orbital periods of months for reasonable satellite 
tidal Q. Inside planet orbital periods of about a week, tidal heating can completely melt the satellite. 
Lastly, we compute an upper limit to the satellite mass loss rate due to thermal evaporation from the 
surface, valid if the satellite's atmosphere is thin and vapor pressure is negligible. Using this upper 
limit, we find that although rocky satellites around hot Jupiters with orbital periods less than a few 
days can be significantly evaporated in their lifetimes, detectable satellites suffer negligible mass loss 
at longer orbital periods. 
Subject headings: 



1. INTRODUCTION 

The high photometric precision of the Hubble Space 
Telescope allows interesting constraints to be placed 
on the existence of massive satellites orbiting tran- 
siting exoplanets. Two systems ha, v e bee n searched 
so far, HP 209458b dBrown et all |200H ) and HD 
189733b (|Pont et al.ll2007[ ). From the transit hghtcurve 
and timing, upper limits of roughly an Earth ra- 
dius and mass have been placed on hypothetical satel- 
lites around these two planets. A number of studies 
have discussed detection t e chniq ues (iBarnes & Fortneyl 
2nnl: jDovle fc Deed IMil iHan fc Han 2002; KioDing 
2009a"b': 'Sarto retti fc Schneideij|199a : .Simon et al..,2007l : 
Szabo et al. 200^ This paper is concerned with the or- 
bits and physical structure of these hypot hetical moons. 

Wi th the launch of the Kepler satellite (jBorucki et al.l 
|2004|) . ~ 10^ stars will be monitored with photometric 
precision sufficient to detect Earth size objects. The es- 
timated hundreds of hot Jupiters (gas giant planets with 
orbital periods less than 1 week) to be found by Kepler 
will provide a large sa mple to inves tigate the existence 
o f Earth-size satellites (ISzabo et al.l [2006). 

iJohnson fc Huggini (|2006( ) note another possible ob- 
servational manifestation of satellites. Interaction be- 
tween a satellite and the planet's magnetosphere may 
give rise to an extended gaseous torus orbiting the planet, 
as in the Jupiter-Io system. Material stripped from the 
satellite's surface and atmosphere would contribute to 
the trans mission spectrum of the p lanet. It has been 
suggested (jVidal-Madjar et al.' 2003) that absorp tion in 
the upper atmosp here of HD 209 458b CBallester et al.l 
2p07; Charbonncau ct al. 2003; iRedfield et all 120081 : 
Vidal-Madjar et al.n2003. . .2004 ) is due to a high rate of 



atmospheric escape which increases the density at large 
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radii. The plasma torus model is an alternative explana- 
tion for such absorption high in the atmosphere. 

There has been interest recently in the properties of 
satellites around exopl anets, mainly focused on the hab- 
itability of moo ns ie.g. iRevnolds et al.ll987t[Sch"ari2006l : 
IWilliams et al.lfT997i ) around planets at large orbital sep- 
aration from the star. Our focus in this paper is on 
planets to be found with Kepler, hence we restrict our 
attention on close-in planets for which the transit proba- 
bility is substantially higher. As we will show, moons in 
such systems will be subjected to intense tidal heating, 
with consequences for the moon's structure, evaporation 
rate, and interaction with its surroundings. 

Two recent studies of th e orbital stability of sa tellites 
are important in this work. iDomingos et al.l (|2006[ ) inves- 
tigated orbital stability for point mass star, planet and 
satellite, for a range of eccentricity of both the planet and 
satellite. They point out that, for small planetary eccen- 
tricity, satellite orbits are stable for as < O.Sajy, where as 
is the semi-major axis of the satellite around the planet 
and an is the planet's Hill radius. As we will show, the 
satellite orbits have large periodic eccentricity induced by 
the star as the stability limit as — 0.5a// is approached. 
Hence, while lo's eccentricity and tidal heating are due to 
interaction with Jupiter and the other Galilean satellites, 
moons of close in exoplanets have a large periodic eccen- 
tricity due to stellar forcing. We show that this periodic 
eccentricity leads to large tidal heating for systems close 
to the parent star. This tidal heating may have a signif- 
icant impact on the structure of the satellite, inducing 
melting of the interior and perhaps reinforcing dynamo 
generation of magnetic fiel ds in the core. 

As pointed out by Bar nes fc O'BrienI ()2002f ). such 
close-in moons, orbiting slowly rotating planets synchro- 
nized to the star, are subject to orbital decay. The tide 
raised on the planet by the moon acts to spin up the 
planet while shrinking the moons orbit, eventually lead- 
ing to the moon impacting the planet. The orbital decay 
rate is sensitive to the tidal Qp of the gas giant planet. 
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iBarnes fc O'Brienl (|2002f) used Qp = 10^ the estimate 
for the Jupiter-Io tidal interaction, leading to the con- 
straint that hot Jupiters cannot have moons larger than 
~ 1O~^M0, where is the Earth's mass. Such tiny 
moons would be undetectable by the transit method. 
However, as we discuss, the appropriate value of Qp may 
be larger by a factor ~ 10^, allowing stable orbits for 
potentially detectable moons with radii > i?^. 

Little attention has been given to the formation sce- 
nario of the massive satellites discussed in this paper. 
ICanuD fc WardI (|2006( ) note that the outer planets of the 
Solar System contain only a fraction 10^'' of the planet's 
mass in the satellite systems. If applied to extrasolar 
planets, this phenomenological scaling implies an upper 
limit for the satellite mass of < M^/SO which would be 
undetectable by the transit method. An alternative pos- 
sibility is for the gas giant to capture a terrestrial planet 
into a bound orbit. 

Satellite orbits in the restricted 3 body problem are 
discussed in section [2] Orbital decay due to the tide 
raised in the planet by the satellite is discussed in section 
[3l Tidal heating of the satellite due to the non-circular 
orbit and its consequences for internal structure are dis- 
cussed in section m Coupled tidal heating and satellite 
orbital evolution are described in section [5l The rate of 
evaporative mass loss from the satellite is discussed in 
section [5] and conclusions are stated in section [T] 

2. SATELLITE ORBITS 

In this section we review the moon's orbit, treating 
moon, planet and star as point masses. This problem 
is referred to as the "main problem in lunar theory" 
([Brouwer &: Clemencd fT96l[ ) . and considerable analytic 
progress can be made. We will give a simple derivation 
of the satellite's orbit, as needed to compute tidal heat- 
ing, for the case where the planet's orbit around the star 
is circular (the "circular restricted 3-body problem" ) and 
the "unperturbed" orbit of the satellite is circular. That 
is, we will compute the small deviations from a circular 
orbit due to the gravity of the star. 

We consider star, planet and satellite of masses ^ 
Mp 3> Ms, respectively. We will ignore the influ- 
ence of the satellite on the star-planet orbit, treating 
it as Keplerian with semi-major axis Op and eccentricity 
Cp. The mean angular motion of the planet's orbit is 
Up ~ (GM^/Op)^/^ = 27r/Pp, where Pp is the planet's or- 
bital period. The satellite's orbit around the planet can 
be treated as nearly Keplerian, with semi-major axis a^, 
eccentricity and mean motion Us — (GMp/a^y^"^ = 
2n/Ps. 

We work in a reference frame with origin at the center 
of the planet and non-rotating axes fixed with respect to 
distant observers. In this reference frame the star only 
contributes a tidal force. Cylindrical coordinates (r^, (pg) 
and (r,, are used for the positions of the satellite and 
star, respectively. The orbits are assumed coplanar for 
simplicity. The equations of motion for the satellite are 
then 
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where £s is the satellite's orbital angular momentum per 
unit mass and TZ = 7?,(rs, c^s, r*, 0*) is the negative of 
the tidal gravitational potential from the star. We can 
simplify the form of TZ using the fa ct that <g; r^, to find 
the l eading order result in r^/r* (jBrouwer fc Clemencd 

[TmH), 
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Given an orbit for the planet-star system, r*(<) and 4>*{t), 
eqUl [2] and [3] can be integrated in time to find the orbit 
of the satellite. 

The parameter space for stable satellite orbits is 
strongly constrained. For close-in planets, the plane- 
tary radius Rp can be a large fraction of the planet's 
Hill radius, an = ap{Mp/3Mt,y^^ , implying that all of 
its bound satellite orbits are signifi cantly perturbed by 
the star's gravity. iDomingos et al.l ([2006 ) found stable 
orbits for as < 0.49a_f/(1.0 - l.Oep - 0.276^). In terms 
of the orbital periods, the stability limit is Pg ^ -Fp/5 
for ep,es <C 1. The dimensionless strength of the dis- 
turbing potential is then TZ/ {GMp/us) ~ {np/ugf' — 
{Ps/PpY < 1/25. Ignoring the cohesive strength, the 
moon must also orbit outside the planet's Roche radius 
at flRochc = 3^/3i?p(Mpi?3/Msi?3)i/3^ where Rp is the ra- 
dius of the planet. For gas giants with mean density pp = 
3Mp/ (AirRp) ~ 1 g cm~'^ and terrestrial satellites with 
mean density ps = 3Ms/(47ri?f) '--^ 5 g cm~^, the Roche 
radius is at or below the planets surface. Hence we can 
consider the satellite disrupted if it hits the planet's sur- 
face. In summary, satellite orbits are bounded at short 
orbital periods by the planet's surface and at long or- 
bital periods by the stability constraint, giving the range 
of possible orbital periods 2'K{Rl/GMpf/^ < Ps < Pp/5, 
for es , Cp <C 1 



The derivation in'Brouwer fc Clemencd (|1961h uses the 
disturbing function formalism and computes the varia- 
tion of the orbital elements Os, Cg, etc rather than rg 
and (ps directly. Here we give a simpler derivation of 
the variation involving Ts and (ps. We consider a cir- 
cular planetary orbit with — ap and p^ — Upt, and 
a satellite orbit described by a circular orbit plus small 
perturbations induced by the star. In detail, we write 



rs{t)=as + 6rs{t) 

(j)s (t) = psQ +nst + 6ps (t) 

£s{t)^alns + 5£s{t). 



(5) 
(6) 

(7) 



Here is a constant phase. Plugging eql5l[6]and[7]into 
eq[Tl[2]and[3l the equations are satisfied at leading order. 
Eq[3] can be immediately integrated to yield 

SUtl ^ 3 /IV y [2(n, - np)t + 2^so] , (8) 
is 4: \ns J 

where we have approximated rip <ti Us to simplify the 
expression. Plugging cq|5]into eq[Tl wc find the following 
forced harmonic oscillator equation for Sr, 

6rs + ii^s^rs Tipfls ( i -f 3 cos [2{ns - np)t + 2Psq\ ) (.9) 
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Fig. 1. — Perturbations to the radius (upper panel) and lon- 
gitude (lower panel) due to the stellar gravitational field. Here 
Sts = Ts{t) — as and 5(f>s = 4's{t) — "n-st- The solid lines are the re- 
sult of numerical integration of eq[T][2]and|3]for Pp/Ps = 20, scaled 
by the factor (Pp/Ps)^. The dashed lines are the analytic results 
including the periodic components from eg llOl and llll as well as the 
time-independent shift to the radius Srs/as = (np/ras)^/2, and the 
resultant change in orbital period. This change in orbital period is 
the cause of the linear beha vior of 5ij> in the lower plot, while the 
oscillatory term from eg 111 I causes the short period variations. 



The first term on the right hand side of eq[S] gives a shift 
in the orbital radius Sts = {l/2)as{np/ns)'^ . This con- 
stant term causes no tidal heating, so we ignore it from 
here on. The second term on the right hand side of eql5] 
causes an oscillatory change in 5rs ■ We can solve for this 
term by writing 5rs{t) = Sr^Q cos ([2(^5 — np)t + 20so]), 
where Srgo is a constant. The cosine factor then cancels 
out of the equation and the amplitude Stso can be found. 
The result is 

Srsit) f ro/ u 

— I > cos [2[ng — Upjt 



24> 



sO\ 



(10) 



Plugging eqlS] and [TO] into eqUJ we find the angle 

2 



11 



sin [2(ns 



,)t + 2(bso]. (11) 



Eg 1 1 01 and 1111 confirm that the forced perturbations vary 
as [Ps/ PpY for small eccentricity orbits. These formulas 
are valid when Ps ^ Pp, and underestimate the pertur- 
bations as the orbit approaches instability. We will use 
eqfTOlandfTTlto compute the tidal heating rate in section 

H 

Figure [1] compares a direct numerical integration of eq. 
[U [2] and [3] against the analytic formulas in eqdO] and [Til 
showing that they agree to the stated accuracy in the 
small parameter Ps/ Pp. 

The discussion so far has examined small perturba- 
tions about circular planet and satellite orbits. The same 
method can be used to find small perturbations, induced 
by the star, about orbits with finite "background" eccen- 
tricities e^o or epQ. A detailed derivation and results can 
be found in Br ouwer fc Clemencd ()1961h . Again orbital 
variations of size {Ps/Pp)"^ are found due to stellar per- 
turbations. Tidal heating of the satellite can arise due 



to either background e^o and Cpo, or perturbations in- 
duced by the star, which exists even for e^o = epo — 0. 
We show in section [Jj that the relative size of these two 
effects scales as Cso or Cpo versus {Pg/PpY- The forced 
variation should provide a lower limit to the tidal heat- 
ing, and may significantly underestimate the heating for 
finite eccentricity orbits at large ap. 

3. ORBITAL DECAY DUE TO TIDES RAISED IN THE GAS 
GIANT PLANET 

Tidal friction is recognized to be an important fac- 
tor in the survival of hypothetical early satellite systems 
for inner So l ar System planets M ercury and Venus (e.g., 
lBurnslll973l : IWard fc Reidlll973D . These studies pointed 
out that tides raised in the planet by the satellite and the 
Sun can cause significant changes in the satellite's orbit. 
The problem is even more severe as the planet is moved 
closer to the star, as the stellar tides are stronger and 
the Hill sphere is closer to the planet. However, a key 
difference occurs if the planet in question is a gas giant 
rather than a rocky body. Gas giants are less dissipa- 
tive than rocky bodies, and hence satellites are relatively 
more immune to destruction by orbital decay into the 
p lanet. 

iBarnes fc O'Brien! (|2002[ ) studied orbital decay of 
satellites around hot Jupiters. They point out that if 
the gas giant planet is synchronized to the parent star, 
then the satellite orbits faster than the planet spins. 
Tides raised in the planet by the satellite then attempt 
to spin up the planet while shrinking the orbit of the 
satellite to conserve angular momentum. The conclu- 
sion of their study was that only satellites with masses 
Ms « could survive being dragged down to the 
planet. Such tiny satellites would be undetectable by 
the transit method. We point out that recent theo- 
retical work on the tidal Qp in gas giants implies that 
orbital decay is fa r less effective than the estimate by 
iBarnes fc O'Brielil ([2002) for the forcing periods of inter- 
est, and hence orbital decay may be ineffective even for 
Earth mass satellites. 

Using the equations from lGoldreich fc Soted (| 1966( 1 . or- 
bital decay over a timescale T will lead to destruction of 
satellites with mass 



M, > M, 



5.4 X lO"-" M„ 



QpPi 



13/3' 



plO/3^ 
^dyn ^ 



(12) 



where Qp is the quality factor for the gas giant planet 
and Pdyn — 2TT{Rp/GMp)^/^ is the dynamical time of 
the planet (2.8 hr for Jupiter). IBarnes fc O'BrienI (|2002[ ) 
evaluated this expression using Qp ~ 10^, the value 
infer red for orbital expansion in the Jupiter-Io system 
( Goldreich fc Soteil il966). as well as circularization of 
the extrasolar planets fe.g.. IWull2003f ). finding 



:5 X lO^^MfT 



Pr, 



Mp 

13/3 



Qp 

105 

2.8 hr 

Pd 



yn 



10/3 



5 Gyr 
T 



^Pp/5 4 days / 

where Mjup is the mass of Jupiter, and we have scaled 
Ps in terms of the maximum orbital period Pp/5. If true, 
this result would imply that Earth mass satellites would 
only be stable outside orbital periods of ^ 40 days, with 
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corresponding smaller transit probability than for the hot 
Jupiters at Pp ~ a few days. 

The physical origin of the tidal Qp in gas giant plan- 
ets has been an outstanding question since at least 
the 1970's. The Qp for Jupiter is constrained to be 
Qp ^ 10^ — 10^ in ord er that lo's orbit expand ed into 
the Laplace resonance (iGoldreich fc Soteri [l966l). Early 
theoretical work by iHubbardI ()1974f ) showed that the 
turbulent viscosity generated by convective eddies re- 
quired for the outward transport of heat would give 
rise to Qp ^ 10^, perhaps explai ning the observed 
value. IGoldreich fc NicholsoiJ (|1977l ) then pointed out 
that turbulent eddies in the planet have long turnover 
times, ioddy, compared to the forcing periods of interest, 
severely decreasing the turbulent viscosity used in Hub- 
bard's calculation. IGoldreich fc NicholsonI (|1977h esti- 
mated Qp ~ 10^'^ for "equilibrium-tide" flow in Jupiter, 
underpredicting the obs erve d tida l dissipation rate by 
a factor of 10^ - 10®. |^ (pOOl revisited Goldreich 
and Nicholson's calculation, revising Qp downward to 
Qp ^ 1012. 

The detailed physics underlying tidal dissipation in gas 
giants may strongly affect predictions of the orbital sta- 
bility of the satell ites stu died in this pa per. Recent cal- 
culations by Ogilvie & Lin (2004) and (|2005D have 
found that when the tidal forcing frequencies are reso- 
nant with inertial waves (waves with frequencies less than 
twice the spin frequency, restored by the Coriolis force) 
that Qp ~ 10^ — 10^ can be obtained, due to the higher 
dissipation rate for short lengthscale inertial waves. Con- 
versely, when the forcing frequency is not in the inertial 
range, no such resonant excitation occurs and Qp will be 
orders of magnitude larger, comparable to the equilib- 
rium tide value Qp ^ 10^^, or perhaps somewhat smaller 
if atmospheric waves can dissipate efficiently. The work 
of Ogilvie and Lin and Wu has demonstrated that res- 
onant excitation of inertial waves is a promising mecha- 
nism to explain the small observed Qp values. A generic 
feature in their results is that the tidal dissipation rate is 
strongly dependent on the forcing frequency. Hence the 
practice of taking a measured Qp from one situation and 
applying it to another may be incorrect if the forcing fre- 
quencies are very different. The consequences of the tidal 
dissipation factor varying with forcing frequency have al- 
ready been explored in the context of solar type stars by 
IQgilvie fc Linr(|200l . 

We now explore the consequences for survival of satel- 
lites of short period exoplanets if the inertial wave hy- 
pothesis is correct. If synchronized, the spin period of the 
gas giant equals its orbital period, Pp. But the orbital 
period of the satellite must satisfy < Pp/5 for orbital 
stability. For a quadrupole tide raised in the planet, the 
forcing period would be shorter than Pp/10, while in- 
ertial waves only exist with periods longer than Pp/2. 
Hence inertial waves cannot be excited in synchronous 
planets by the satellite and Qp ~ 10*^ does not apply. 

The existence of exomoons around close-in gas giant 
exoplanets with Pp a couple months provides a test 
of tidal dissipation theory. From eqll31 in order for a 
satellite of mass Ms to survive for a time T requires 



h > gp,crit = 2 X 10^ 



M 



Jup 



Mr, 



Pp 4 days 

5Ps Pr> 



13/3 



P 



dyn 



2.8 hr 



10/3 



T 



5 Gyr 



(14) 



This lower limit on Qp is well above the canonical value 
Qp ~ 10^ — 10^, and thus could confirm the expected 
large difference in Qp between forcing in and out of the 
frequency range of inertial waves. Outside the inertial 
range, the equilibrium tide calculation gives Qp ^ 10^^ 
(IGoldreich fc NicholsonI fl977l : Iwjll2005h . so that Earth- 
size satellites of hot Jupiters would be immune to de- 
struction by orbital decay over Gyr timescales. Note 
however, that even if Qp ~ 10^ — 10^, Earth size satel- 
lites could exist for longer period orbits Pp ^ a couple 
months. 

4. HEATING OF THE SATELLITES OF CLOSE-IN 
EXOPLANETS 

This far we have studied two facets of satellite orbits. 
First the variations induced in the satellite orbit by stel- 
lar gravity were reviewed. Then we showed that satellite 
orbits may be far less susceptible to orbital decay than 
previously thought if Qp ^ 10^ — 10^ for forcing outside 
the inertial frequency range, as occurs for a satellite or- 
bit around a synchronized planet. Given that long-lived 
satellites may exist around even short-period exoplanets, 
we now investigate the consequences of their orbital vari- 
ation on tidal heating. To motivate our study, we first 
review tidal heating in lo. 

Based on the significant forced eccentricity of lo, 
iPeale et all (|1979D predicted widespread surface volcan- 
ism, wh ich was soon confirm ed by images from Voy- 
ager 1 (jMasurskv et al.l Il979l ). The total dissipation 
rate for tidal for cing of a homogeneous, incompressible , 
elastic sphere is (jPeale fc Cassenlll978t iPeale et aLlll979l: 
lWisdomll200l 



E- 



42 TTpin^R'^ejo 

19 ^J.sQs 



(15) 



where /i^. ~ 6.5 x 10^^ dyne cm"^ is the rigidity, and 
Qs is the tidal dissipation coefficient of lo. For lo, 
Ps — 1.7days, ps — 3.5g cm^^, and Rs = 1800 km. 
The eccentricity o f lo was determined to be Cso — 0.0043 
(jPeale et al.lll979D . resulting in a dissipation rate for lo 
of 



Eio — ■ 



1.6 X 10^1 



Qs 



erg s 



(16) 



IPeale et all ()1979D used a fiducial value Qs — 100, based 
on laboratory studies on rock samples. For Qs — 100, 
eq ll6l gives a heating rate several ti mes the r adioactive 
heating rate estimate for the Moon fPeale et al. 1979). 
Tidal heating is widely recognized as an important effect 
in the thermal structure of lo. 

To assess the role of tidal heating for the satellites 
around hot Jupiters, we first ask what would happen to 
an lo-like satellite orbiting around a Jupiter- like planet, 
which itself is in 4 day orbit around a solar-type star. 
From eq llOl and 111! the perturbations correspond to 
eccentricity variations ^ {Ps/PpY ^ 1/25 — 0.04. 
For orbits near the stability limit, the forced eccentric- 
ity is larger than the eccentricity of lo by a factor of 
0.04/0.0043 = 9. The orbital period of the satelhte must 
be smaller than 4 days/5 = 19hr, shorter than lo's or- 
bital period by a factor of ~ 2.1. Plugging into eq[T5l the 
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higher eccentricity and shorter orbital period combine 
to increase the heating rate by a factor ^ W^. Clearly 
tidal heating will be important for satellites around hot 
Jupiters even for satellite size much smaller than lo. 

Tidal heating rates for Earth-like satellites are sen- 
sitively dependent on the tidal Qs- The value 
QEarth-Moon ~ 10 for the Earth, as derived from Earth- 
Moon orbital evol ution, is likely due to dissipation in 
shallow seas (e.g. lJeffrevslll970l) . Heat from such dis- 
sipation depends on the uncertain ocean depth and to- 
pography, and is easily radiated away. We focus instead 
on heating of the interior, especially for early stages to 
determine if the satellite is melted. Decay of seismic 
body waves with periods of 10^ — 10"^ s in the Earth finds 
Qsarth seismic ~ 300 - 500 (e.g. iRoult ct al.l l200l . but 
extra polation of this data to ~ 10^ — 10^ s is uncertain 
(e.g. iBurndllQT?! : lEfroimskv fc Lazarianll2000D . For the 
solid phase, we will use a fiducial value of Qs = 100, 
assuming moderately increased dissipation when the pe- 
riod is extrapolated over 10 — 100 longer timescales than 
is measured. As we will show, the large tidal heating 
rates for a solid satellite with Qs ^ 100 imply that satel- 
lite is likely melted on a short timescale. The subsequent 
tidal heating rate depends on the appropriate Qs for the 
liquid state, which is uncertain. For a completely molten 
satellite, Qs may be much larger. For instance, the Pre- 
liminary Reference Earth Model, constructed from seis- 
molo gical data, uses Qs ^ 10^ for th e Earth's liquid outer 
core (jDziewonski fc Andersonlll981[ ). 

The tidal heating rate in eg fT51 derived bv'Peal e et al.l 
([1979) is for the case of a Keplerian satellite orbit with 
eccentricity e^o. This is only strictly valid for a satel- 
lite orbit unperturbed by external influence. Here we 
are concerned with satellite orbits which may be signif- 
icantly perturbed by stellar gravity. In other words, we 
require a tidal heating rate for satellite orbits given by 
the solution to the restricted 3-body problem (for small 
satellite mass), rather than for a Keplerian orbit. The 
orbital variation due to stellar gravity was summarized 
in section [21 

To discuss the tidal heating rate for satellites perturbed 
by stellar gravity we will derive an analytic est imate by 

S iluggi ng eg fTUl and [TT] into the formula found in I Wisdom! 
200l : 



E = 



PshsRs 



dn U 



dU' 



(17) 



gs J dt 

This more general formula allows an arbitrary orbit, and 
is not specialized to the case of an unperturbed Ke- 
plerian orbit. Here hg = (5/2)/(l -I- Id^isf^PsgsRs), 
gs — GMs/R^, U is the tidal potential in the satellite 
due to the planet, U' is the tidal potential including the 
dissipative lag, and the integral extends over the surface 
of the satellite. We perform several operations to sim- 
plify eq[T71 First, we expand U in spherical harmonics 
and perform the angular integration. Next, we expand 
to leading order in the lag time (tIsQ^)"^. We assume 
the satellite is synchronized when evaluating the libra- 
tion term. Lastly, we plug in the small deviations from a 
circular orbit found in eg fTUl and [TT] We find the result 

■ _28n GMlhsRlns ( npV 

Qsal \-s) ■ ('^^ 
Note the difference between eq.s [151 and [181 Eq[T5l is 
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Fig. 2. — Satellite tida l heating rates (solid lines) as a function 
of orbital period. Eg 1191 was used for the heating rate with ps = 
5 g cm~'^, Qa = 10 and Ps = Pp/5. Insolation rates (dashed 
lines) for a Sun-like star are shown for comparison. 



proportional to e^g, which is assumed constant for the 
unperturbed Kepler orbit, while eqlTSlis proportional to 
{up/us)'^, the rms eccentricity squared induced in the 
orbit by the stellar gravity. No terms involving Cso 
appear in eg |18l as we have perturbed around a Kep- 
lerian orbit with e^o = 0. Alternatively, a more in- 
volved calculation would be to perturb around an orbit 
with finite e^o, in which case presumably the terms in 
eo ll5l and eg ll8l would appear summed together. That 
is, plugging the Keplerian result Sr — — CsoOs cos{nst) 
and 6(f)s = 2eso sm{nst) into eg |17l would yield eg |15l 
Equating these two formulas we find a critical eccentric- 
ity es,crit — l-6{Ps/Pp)^ ^ 0.06 above which the term 
in egllSI would dominate and below which the term in 
eq ll8l would dominate. The eccentricity e^o would de- 
pend on the initial conditions as well as subsequent tidal 
dissipation in the satellite and planet which would act to 
decrease Csq. Henceforth we ignore e^o and use the tidal 
heating rate in eqll8l Since we have ignored possible fi- 
nite Cso, we may underestimate the heating rate if eso is 
large. 

To numerically evaluate eg ll8[ we again use fis — 6.5 x 
10^^ dyne cm^^ for the rigidity, and scale all quantities 
to Earth values. Since the heating rate increases as Ps 
decreases, we scale expressions to the minimum heating 
rate using Pg = Pp/?>, finding 



^solid 



3.2 X 10^^ erg s'^ / p. 



Qs V 5 g cm 3 

RsV daysV /Pp/5 
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(19) 



for a solid satellite and 

2.7 X 10^^ erg s^^ 



-E'liquid — 



Equivalcntly one could use eg 1291 which is simpler in practice. 
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Fig. 3. — Effective temperature due to tidal heating (solid 
lines) or insolation (d ashed line). Solid lines show temperature 
determined from eg 12 II while dashed lines represent the equilibrium 
temperature due to insolation, and reradiation over 4tt steradians. 
Tidal heating determines the temperature when the solid line is 
above the dashed line. 



4 daysV fPp/5 



Pr, 



Ps 



(20) 



for a satellite which has been completely melted {^s — 
0). Eqliniis compared to insolation in figlH While the 
prefactors in eqlTH] and [201 are comparable, their scaling 
with Rs differs, and the tidal Qs may be quite different 
for molten and solid satellites. If the tidal Qs in the liquid 
phase is far larger than for the solid phase, then the rate 
of tidal heating could drop drastically upon melting. 

Would Earth-size satellites around hot Jupiters be 
melted by tidal heating? Eqlin] implies a heating rate 

E^oiiA - 10^* erg s- 
latent he at L = 10^° erg ap i 

pressure (jStevenson et al.iri98 3^. melting would occur in 
a few hundred years. The tidal heating rate, per gram, 
scales as so that heating is more intense for larger 
satellites. Assuming the radioactive heating rate of the 
Earth, ~ 10^° erg s^^, tidal heating dominates radioac- 
tive heating for Pp < 150 days. Radioactive heating 
will eventually dominate for satellites far smaller than 
an Earth radius. The rmal models of the young Earth 
(jStevenson et al.l [19831 ) typically show a completely liq- 
uid core for E > 10^° erg s~^; the tidal heating rate is 
greater by many orders of magnitude. We conclude that 
the Earth-size satellites of hot Jupiters, observable by 
their transits, would be largely molten, except perhaps 
for a thin layer ne ar the surface, cooled by radiation (e.g., 
iPeale et al.l[l979l ). HQs drops precipitously for com- 
plete melting, the temperature could be regulated near 
the melting point. 

The effective temperature for the tidal heating energy 
flux may be found by equating the heating rate in eq ll9l 
to the blackbody coohng rate AnR^aT^ff, giving 



^ for a 4 day orbital period. For a 
erg ap propriate for iron at high 



Teff = 1800 K 
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Tidal heating dominates the temperature structure, even 
at the surface, when the temperature in eq l21l is larger 
than the equilibrium temperature implied by insolation. 
Figure [31 compares these temperatures for Qs — 10^, for 
a zero-albedo surface with uniform temperature over the 
entire surface. Clearly for Qs = 10^, Earth-sized objects 
are raised above the melting temperature at zero pressure 
(Imeit ^ 2000 K) for either iron or rock compositions 
inside a critical orbital period, which we estimate as 

/ 2000 K y/' f Ps 



^'p,mcit = 3.7 days 



5 g cm 



Rs 

R(f 



(22) 



Qs 

Hence, at the orbital periods of the hot Jupiters (a 
few days), satellites Earth size or larger are completely 
molten for Qs = 10^. 

Vigorous tidal heating has consequences for the satel- 
lite magnetic field. Fluid motions in the conducting core 
driven by tidal heating may generate magneti c fields 
through dynamo action. jStevenson et al.l ()1983l ) derive 
a magnetic field for the Earth by equating the energy 
available for dynamo generation to the Ohmic dissipa- 
tion rate of the field, which yields the scaling cx E, 
where E is the tidal heating rate. As the dissipation rate 
in eq[20| is larger than the Earth's cooling rate by a fac- 
tor ^ 10^ /Qs, our results imply that fields larger than 
Earth's can be created if Qs < 10^. Intrinsic magnetic 
field can decrease the mass loss rate of the atmosphere. 
Magnetic interaction between the satellite and planetary 
magnetosphere may be of the Jupiter-Ganymede type 
(jKivelson et al.|[l997[ ) in which the standoff distance is 
determined mainly by magnetic stresses. 

Stu dies of habitability o f satellites around gas giants 
(e^ [Reynolds et al.l[l987t iScharJ [20061 : [Williams et aH 
[1997f ) have invoked tidal heating to maintain plate tec- 
tonics, create subsurface oceans, etc. They computed 
satellite tidal heating for finite eccentricity, assuming it 
to be a free parameter, or pumped to large values by in- 
teraction with another satellite. The mechanism for tidal 
heating in this paper, periodic eccentricity forced by the 
star, is likely unimportant for satellite habitability, since 
it is only large for systems very close to the star, where 
water is already in liquid or vapor form. Tidal heating is 
of increasing importance, relative to insolation, for host 
stars of smaller mass since the stellar luminosity drops 
rapidly. 

5. ENERGETICS 

Tidal dissipation in the satellite takes energy out of 
the satellite orbit and deposits it in the satellite body 
as heat. In this section we show that the energy lost 
from the satellite orbit by heating is exactly balanced by 
energy gain from work by the stellar gravity. The tidal 
friction acts to create a lag in the satellite orbital veloc- 
ity. This lag has the correct phase to allow energy input 
to the satellite's orbit at a rate sufficient to power the 
tidal dissipation in the satellite. Ultimately, the energy 
reservoir for tidal heating of the satellite must then be 
the star-planet orbit. 
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Tidal friction acts to exert a velocity dependent accel- 
erati on on the reduced ma ss of the satellite-planet sys- 
tem (iMardling fc Linll2002l . eq.4) with radial and angular 
components given by 



fTF) _ 



-70 rs 
-70 rs 



(23) 
(24) 



Here 70 = {9/5){nshs/Qs)imp/ms){Rs/as)^ is the fric- 
tional drag coefficient, which is much smaller than the 
orbital mean motion ris . Again we choose a circular back- 
ground orbit plus small perturbations as in eqlSl [S] and 
H EqESl and [M] become 

.(TF) 



(TF) 



-70 Sf-s 

= -70 as (us - fls + 64>< 



(25) 
(26) 



Assuming that the satellite spin has already reached the 
synchronous state, we may set ^Ig = Ug in eg 1261 

Adding the tidal friction acceleration in eqUH] and [221 
into eqI5] and [51 and again ignoring the constant radial 
force term, we find equations which have the form of a 
damped, driven oscillator. In the limit 70 ^ Ug we find 
the solutions 



Qg 



= 2 



sin [2{ns - np)t + 20so] 



35 ^ , 

-70 ^-j cos[2(n. 



np)t + 2(j)go] (27) 



cos [2{ng — np)t -t- 2(j}gQ\ 



173 / \ . 

— 70 -) sm[2(n. 



np)t + 2(j)go]. (28) 



The first terms in eq[27| and [^Sj agree with eq[Tni and [TTl 

The second terms proportional to 70 describe a lag in the 
satellite orbit due to tidal friction. 

The energy lost from the orbit is found by computing 
the work done by the tidal friction force. Since the tidal 
friction force already contains the small parameter 70 ^ 
Ug, we can ignore terms of order 70 in eq |27l and l28l For 
reduced mass jigp = MgMp/ [Mg -f Mp) ~ Ms, and time 
averaging so that cos^,sin^ ^ 1/2, we find 



E. 
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70 iSr-g 



313 
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This is just a rederivation of eqlTSl 

The increase in satellite orbital energy due to work 
done by the stellar gravitational force in eq[5| and [5) is 



f-^sp 



-nlag {Sf-g cos [2{ng - np)t + 2(/)so] 



s54>g sin [2{ng - np)t + 20^0. 



(30) 



The terms independent of 70 in eq |27l and [2S1 time aver- 
age to zero when inserted in ea l30l Inserting the terms 
proportional to 70 into eq[30]and time averaging we find 

/ o \ 2 



(31) 



exactly equal and opposite to the energy lost to tidal 
heating in eql29l 

We have found that energy lost from the orbit is replen- 
ished by work done on the orbit by the stellar tidal force. 
Hence the satellite orbit is stable on long timescales to 
tidal effects arising from dissipation in the satellite. The 
energy reservoir powering the tidal dissipation in the 
satellite must then come from the star-planet orbit. For 
star-planet orbital energy 
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GM^Mp 
2ar, 



-1.8 X 10^" erg 
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(32) 



and dissipation in a liquid satellite (eg (20)) . the star- 
planet orbit is immune to orbital decay on a timescale 
T if 



Q.>480 ( ^ 



-2/3 



4 days 



13/3 



( Mp 
\Mjup 
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Rs 

Rcf 
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10 
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(33) 



This critical value of Qg for orbital decay is comparable 
to the value expected for a solid Earth-like satellite (see 
section [4]). In the more likely case that the satellite has 
been melted by tidal heating and Qg ^ 10^, the star- 
planet orbit is immune to decay. 

6. EVAPORATIVE MASS LOSS 

The large tidal heating rates in eg fT^ and large in- 
solation from the parent star, and low escape speed raise 
the possibility of enhanced rates of mass loss from the 
satellites studied in this paper. The rates of mass loss 
are difficult to predict, and may depend on many factors 
such as stellar wind or magnetospheric erosion of the at- 
mosphere, and the flux of vapor from the surface to the 
exosphere. Here we discuss one possible scenario, of ther- 
mal evaporation from the solid surface of a satellite with 
a negligible atmosphere (such as the Moon o r Merc ury) . 

iCameronI H1985D and iFeglev fc CameronI (|1987[ ) dis- 
cussed evaporative mass loss from Mercury in the proto- 
planetary nebula as a means of understanding Mercury's 
high mean density. The high temperatures in the nebula 
allow enhanced rates of evaporation of low(er) density 
silicates from the solid surface into gaseous form. If this 
vapor is removed from the atmosphere, the mass fraction 
of iron in the body, and hence the mean density, will in- 
crease. Due to the exponential temperature dependence, 
this process is only effective at high temperatures. We 
find that for surface temperatures due to tidal heating, or 
the equilibrium temperature due to insolation, that sig- 
nificant ma ss loss can occu r. We will follow the simple 
estimate in ICameronI (|1985D to estimate the importance 
of this process. 

The number fiux for thermal evaporation from a sur- 
face, in cm~^ s~^, is given by 



^cq(r) 

y/2'KmkhT 



(34) 



where Pcq{T) = PqIO is the equilibrium vapor 

pressure above the surface, m is the mean molecular 
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weight (in grams) of the products, T is the teraper- 
ature, and Pq a nd Tp describe the equihbrium pres- 
sure. ICameroiil ()1985f ) considers an illustrative exam- 
ple of MgSiOafliq) = Mgf gas) + SiO(gas) + 02(gas) (see 
iFeglev fc Cameronl (jl987) for more detailed calculation 
of silicate magmas of chondritic composition) , for which 
Po = 10^3.176 (jyj^g cm-'^, To = 24605 K, and m ^ 
5.5 X 10~^^ g. If the surface temperature is near the 
equilibrium temperature for a zero albedo surface (ne- 
glecting tidal heating), radiating with the same efficiency 
over An the equilibrium temperature versus orbital ra- 
dius for a solar type star is T = 5777 K(i?0/2ap)^/^ = 

1300 K(lOi?0/ap)i/^ As seen in figure H this temper- 
ature including only insolation, and ignoring tidal heat- 
ing, may underestimate the satellite surface temperature 
close to the star. Approximating the exponential as a 
power law (valid for small changes in Up), we find a mass 
flux 



m(?!)~1.3 X 10^" g cm^2 g-i 



10i?f. 



21.5 



(35) 



from the surface. If the vapor is lost from the atmo- 
sphere, this translates into a mass loss rate 



M^ = 3.3 X lO^-^M® Gyr^i 



10i?c 



21.5 



Rs 

Ra 



(36) 



In 5 Gyr, an Earth mass can be lost for Up < 7ARq, 
where T > 1500 K. Hence satellites originally of Earth- 
size, and hence detectable by their transit signature with 
Kepler, may become undetectable due to erosion. Eq[5Sl 
implies the rate of decrease of the radius is independent 
of radius. Erosion does not slow for small satellites. 

Our estimates give an upper limit to the erosion rate, 
and show that large mass loss can occur over the lifetime 
of the satellite. The exponential dependence on surface 
temperature suggests a critical orbital radius (of planet 
around the star) inside of which satellites can be sig- 
nificantly eroded. Our conservative estimate finds the 
critical orbital period to be around 2.4 days (6.6 days 
using the day-side equilibrium temperature). 

7. CONCLUSIONS 

* We note that if the atmosphere is thin (surface pressure <C 
1 bar), the day-side temperature may be larger by a factor of 2^/^. 
This will cause the critical ap for large mass loss to increase by a 
factor of 2, and the critical planet orbital period to increase by a 
factor 2^/^. Hence a critical orbital period inside of which large 



In this paper we have studied the orbits, tidal heating 
and evaporative mass loss rates from Earth-sized satel- 
lites orbiting hot Jupiters. After reviewing the pertur- 
bation to satellite orbits due to the stellar gravity, we 
derived four main results. First, we have shown that or- 
bital decay due to tides raised in the planet by the satel- 
lite may be much less efficient that previously thought 
since gas giants are likely very weakly dissipative at the 
forcing frequencies of interest. Consequently, even Earth- 
size satellites may be stable around hot Jupiters over 
Gyr timescales. Second, large tidal dissipation rates are 
induced in the satellite due to its forced orbital varia- 
tions, likely melting all but perhaps a thin surface layer. 
Third, we found that the satellite orbit does not evolve 
secularly due to tidal dissipation in the satellite, as the 
stellar gravity does work at a rate to keep the orbital en- 
ergy constant. Lastly, the estimated upper limit to mass 
loss, valid for sufficiently thin atmospheres, indicates a 
critical orbital period around a Sun-like star (2.4 days 
for a thick atmosphere, 6.6 days for a thin atmosphere) 
for Earth-like satellites orbiting a hot Jupiter. Inside 
this critical period significant erosion can occur, reduc- 
ing satellites to a size undetectable by upcoming transit 
observations with the Kepler satellite. Outside the crit- 
ical period, evaporative mass loss will become negligible 
so that detectable Earth size planets can survive. 
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mass loss occurs depends sensitively on the 2^^^, since it occurs in 
an exponent. Our assumption of equal day-night temperatures is 
conservative. 
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